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'. ABSTRACT 



We have used composite spectra generated from more than 22000 QSOs observed 
in the course of the 2dF and 6dF QSO Redshift Surveys to investigate the relationship 
between the velocity width of emission lines and QSO luminosity. We find that the 
, velocity width of the broad emission lines H/3, H7, Mg II, C III] and C IV are correlated 

■ with the continuum luminosity, with a significance of more than 99 per cent. Of the 

O ! major narrow emission fines ([0 III] A5GG7, [O II] A3727, [Ne III] A3870 and [Ne v] 

CO . A3426) only [O III] exhibits a significant correlation between line width and luminosity. 

' Assuming that the gas is moving in Keplerian orbits and that the radius of the broad 

(-H I line region is related to the QSO continuum luminosity, we use the velocity widths 

Q^. of the broad lines to derive average black hole masses for the QSOs contributing to 

I ' the composite spectra. The resultant QSO mass-luminosity relationship is consistent 

O I with M oc i^O-97±o.i6_ -y^g gj^j that the correlation between line width and redshift, 

if present, must be weak, and only C IV shows significant evidence of evolution. This 
enables us to constrain the redshift evolution of the black hole mass-luminosity ratio to 
be ~ (1-l-z)'' with /3<1, much less than the ~ evolution seen in QSO luminosity 

evolution. Assuming that the motion of the broad line region gas is Keplerian and that 
its radius depends on the QSO luminosity, our models indicate that the observed weak 
redshift dependence is too small for the observed QSO luminosity function to be due 



' to the evolution of a single long-lived population of sources. 

Key words: galaxies: active ~ quasars: general - quasars: emission lines - galaxies: 
stellar content 



1 INTRODUCTION 

The vast energy budget required by active galactic nuclei 
(AGN) is probably provided by the accretion of matter onto 
a supermassive black hole (BH). Evidence that such objects 
can be found in AGN is provided by the Keplerian motions 
of the mega-masers observed in some sources (e.g. Miyoshi 
et al. 1995), the broad FeKa emission line observed in the X- 
ray spectrum of some sources (e.g. MCG-6-30-15; see Tanaka 
et al. 1995; Wilms et al. 2001) and the reverberation map- 
ping experiments (e.g. Blanford & McKee 1982; Netzer & 
Peterson 1997; Wandel, Peterson & Malkan 1999; Kaspi et 
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al. 2000, hereafter KOO; Peterson et al. 2000). The gas in the 
broad line region (BLR) resides close to the central contin- 
uum source and can therefore be used as a probe of the cen- 
tral mass. In particular, if the BLR is gravitationally bound 
and in near-Keplerian orbits, it should be possible to esti- 
mate the central mass, M, from the mean radius, r, and the 
velocity dispersion, v, of the emitting gas. 

Evidence that the velocity dispersions of several emis- 
sion lines, in particular H/3, are indeed dominated by Keple- 
rian motions is presented by Peterson & Wandel (1999) and 
Onken & Peterson (2002). Velocity widths measured either 
from the rms spectra or the mean spectrum, combined with 
reverberation mapping measurements of the BLR radius, 
are now available for more than 34 sources (KOO, and refer- 
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ences therein; Onkcn & Peterson 2002; Woo & Urry 2002, 
and references therein). Although there are some problems 
with this method (see Krolick 2001 for a detailed discus- 
sion), comparison of masses obtained in this way with those 
obtained from the bulge stellar velocity dispersion have been 
found to be in good agreement (Ferrarese et al. 2001, Geb- 
hardt et al. 2000, Woo & Urry 2002). 

Reverberation mapping is extremely time consuming, 
requiring that source light curves are well sampled over a 
number of years and that the sources themselves are vari- 
able. As a result of these constraints, reverberation map- 
ping has only been used to estimate the virial mass in a 
small number (<40) of relatively low redshift objects. It is 
not practical to expand these studies to obtain masses for a 
larger number of sources spanning a larger range in redshift 
and luminosity. 

The available reverberation mapping results show a 
clear correlation between the radius of the H/3 emitting 
region, rap, and the monochromatic luminosity at 5100A 
such that rH/3 oc AL^iqq (Kaspi et al. 2000; Peterson et al. 
2000; McLure & Jarvis 2002; Vestergaard 2002; Netzer 2003) 
where 0.5 < 7 < 0.7. The uncertainty on the slope 7 reflects 
the large scatter in the data as well as the statistical method 
used for the linear regression. These issues have been dis- 
cussed, in detail, by Vestergaard (2002), McLure and Jarvis 
(2002), Maoz (2002) and others. The most recent paper on 
this topic, by Netzer (2003), also addresses these issues and 
comments on the use of the so called "photoionization ra- 
dius" (see Wandel et al. 1999). Regardless of the exact value 
of 7, such a relationship can be used to estimate rn/s in every 
source of known luminosity and to obtain the BH mass by 
combining this size with the observed H/3 line width. This 
method has been referred to as the single epoch mass esti- 
mate. 

Reverberation mapping experiments have concentrated 
on the H/3 emission line. This has restricted the use of the 
"photionization method" to the measurement of BH mass in 
relatively low redshifts {z <0.9) and low luminosity objects, 
although there are some cases of high-z quasars where the 
H/3 line has been observed in the infrared (Shields et al. 
2002 and references therein). The obvious extension to high 
redshift sources is to use the UV luminosity as an estimate 
of r and the width of a certain UV line (e.g. Mg II, C IV 
or Lya) as a measure of the gas velocity. Unlike the case of 
H/3, where the location of the line emission (rna) is known, 
there is no direct measure of the location of the UV lines 
except in a handful of sources. Thus, the only solution is to 
calibrate the known relationship for H/3 against other lines 
assuming the two must result in the same derived BH mass. 
This must be done for sources where rna is known and where 
the velocity width of the UV line and the UV contirmum 
luminosity can be measured. Such intercalibration has been 
performed, successfully, for C IV by Vestergaard (2002) and 
for Mg II by McLure & Jarvis (2002). Thus, both C IV and 
Mg II can be used to obtain single-epoch BH masses. 

Although the results from both Vestergaard (2002) and 
Mclurc & Jarvis (2002) are encouraging, their correlations 
display a large amount of scatter. This may well be intrin- 
sic, rather than due to poor signal to noise. One possibility 
(Young et al. 1999; McLure & Dunlop 2002; Smith et al. 
2002) is that some (or even most) BLRs display a flattened 
disc-like geometry. In this case, the emission line widths 



would be a function of the viewing angle, with larger widths 
observed at larger inclinations, resulting in the velocity dis- 
persion (and hence mass) of a face-on system being underes- 
timated. In addition, the KOO original relationship which is 
the basis of all further work, is hmited to ALa(51001) < 10*'' 
ergs s~^. Thus, mass estimates for the highest luminosity 
sources are uncertain because of the necessary extrapolar- 
tion. 

In this paper, we make use of high signal-to-noise ratio 
(SNR) composite spectra from the 2dF QSO Redshift Sur- 
vey 2QZ (2QZ; see Groom et al. 2002) to test the velocity- 
luminosity relationship in QSOs and to obtain new BH mass 
estimates. The use of composite spectra minimizes the scat- 
ter in the relation between luminosity and velocity width 
caused by effects such as inclination angle. It also minimizes 
the influence of objects with extremely small or extremely 
large line widths. The broad range in luminosity and red- 
shift allows us to investigate correlations of velocity width 
and BH mass with redshift and luminosity. 

In Section 2 we briefly describe the primary spectral 
database and outline the procedure used to obtain the com- 
posite spectra. We discuss the method used to measure the 
emission line widths in Section 3 and derive relations be- 
tween velocity width and luminosity for a number of emis- 
sion lines in Section 4. In Section 5 we extend the work to 
estimate black hole masses from the UV broad emission line 
widths and the QSO luminosity, calibrating them against 
the reverberation mapping results for H/3, and discuss the 
implications of these results for both QSO evolution and 
BLR structure. Our conclusions are summarized in Section 
6. 



2 DATA 

2.1 Observations 

The data used in our analysis is taken from the 2dF and 
6dF QSO Redshift Surveys (2QZ, Groom et al. 2001; 6QZ 
Groom et al. 2003, in preparation). QSO candidates were 
selected from UK Schmidt photographic material over the 
range 16.0 < 6,j < 20.85 using their ubjr colours. Photo- 
metric errors in the bj band were typically 0.1 mag. Objects 
brighter than bj = 18.25 were observed with the 6dF facil- 
ity on the UK Schmidt over the range 3900 7600A with a 
spectral resolution of llA. Those fainter than 6j = 18.25 
were observed with 2dF on the AAT with a spectral range 
of 3700-7900 A and a spectral resolution of 9A. Typical red- 
shift errors are <Tz = 0.003. The median SNR was 15 for the 
6dF spectra and 5 for the 2dF spectra. In total 22041 inde- 
pendent QSO spectra from the 2dF and 6dF observations 
were included in the dataset. 

Absolute magnitudes (Mb) were computed from the 
observed photographic 6j magnitude, after correction for 
Galactic extinction (Schlegel et al. 1998), using the K- 
corrections found by Gristiani & Vio (1990). Throughout 
we assume a flat cosmological world model with = 0.3, 
Ao = 0.7 and Ho = 70 kms"^ Mpc"^ 
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Figure 1. The QSO composite spectrum composed of all the QSOs used in our analysis. Emission lines labeled with solid lines were 
included in our study, lines labeled with dotted lines were excluded. 



2.2 Composite spectra 

Pull details of the creation of the composite spectra and line 
profile fitting is given by Groom et al. (2002). Briefly, com- 
posite QSO spectra (excluding BAL QSOs) were created in 
discrete absolute magnitude (AMb = 0.5 mag) and redshift 
[Az = 0.25) bins. A further series of composite spectra were 
obtained from the QSOs binned in absolute magnitude only, 
as well as a composite contained all QSOs. The numbers of 
QSOs contributing to each composite axe given in Table 1 
in Groom et al. (2002), typically ranging from 100 to 1000 
objects. 

The individual spectra were not flux calibrated and each 
individual spectrum was normalised by a continuum fitted 
between known emission features. Due to the limited wave- 
length coverage of these spectra it was necessary to treat the 
large Fe II features underlying the Mg II A2798 emission line 
as continuum. After continuum normalization the spectra 
were shifted to the rest frame, interpolating linearly onto a 
uniform scale of lAbin"'^. Finally the composite spectrum 
in each Mb , z bin was produced by taking the median value 
in each wavelength element. The median z and Mb of the 
contributing QSOs was also determined for each wavelength 
element 

The composite formed from all 22041 QSOs is shown in 
Figure 1. A number of emission features can be clearly iden- 
tifled, including 3 narrow (forbidden) lines, 9 broad (per- 
mitted) emission lines and one semi-forbidden line (C III] 
A1909). The features identified in Figure 1 with solid lines 
are suitable for detailed study because they exhibit large 
equivalent widths (e.g. Lya, C IV A1909), are relatively free 
from contamination by other emission lines and are visible 
over a large redshift range. Ha or Lyd (Fig.l; dotted lines) 
were not included in our study as, due to the redshift distri- 
bution of the sample, too few spectra contain these emission 
lines for meaningful conclusions to be drawn. 



2.3 Limitations 

Although the use of composite spectra provides a good pic- 
ture of the ensemble average properities of the QSO popula- 
tion, there may be small subsamples of the population e.g. 
weak-lined QSOs for which the any general trends identified 
below do not apply. 

The use of composite spectra may also hide magnitude- 
dependent biases against QSOs with, in particular, weak 
broad emission lines which would prove more difficult to 
identify in the original low SNR spectra. Fortunately, al- 
though the spectroscopic identification of QSOs in the 2QZ 
is a function of apparent magnitude, the identification rate 
is generally high. Groom et al. (2003) provide an empirical 
fit to the overall spectroscopic identification rate (/s) as a 
function of apparent 6j magnitude: 

exp(&j - 20.39) 
(l-0.89)-o-92 

For the 6QZ survey, the spectroscopic identification rate is 
more-or-less constant at /s = 0.97. Thus the identification 
rate for the survey remains above 90 per cent for all but the 
faintest 0.75 mag range of the combined sample. ^ Even for 
the faintest magnitude, the identification rate does not go 
below 80 per cent. 

Finally, we note that as redshift increases, the observed 
spectrum corresponds to a different portion of the rest- frame 
spectrum. The continuum normalization of the individual 
spectra that make up the composite will therefore be defined 
by different regions of rest-wavelength. This has been largely 
circumvented by fixing continuum bands in rest-wavelength 
near the lines of interest and interpolating between them to 

t Based on repeat observations of sources for which we initially 
failed to secure a spectroscopic identification, we also confirmed 
that the fraction of QSOs amongst the unidentified sources is 
similar to the fraction of QSOs amongst the identified sources. 
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define the continuum beneath emission lines. Comparison 
of composite spectra generated for different redshift bins 
within the same luminosity range do not reveal any obvi- 
ous evidence of such systematic redshift effects. It is pos- 
sible, however, that low level systematic redshift biases arc 
present in our data, but we do not believe that they would 
be sufficiently strong to alter our main conclusions. 



3 EMISSION LINE WIDTHS 
3.1 Fitting 

The majority of the strong emission lines in the composite 
spectra used in this study are blended with other weaker 
emission lines, usually from different elements. Additionally, 
permitted emission lines such as H/3 and H7 often exhibit 
both a broad and a narrow component which arc emitted 
from physically distinct regions and we therefore attempted 
to separate these components. 

To measure line widths, the local "pseudo-continuum" 
on either side of each feature was fitted with a straight line 
and subtracted. We modeled the overlapping lines contribut- 
ing to each spectral feature using multi-component Gaussian 
fits. We note, however, that assuming a Gaussian form for 
the features in our spectra is a gross over-simplification, and 
thus the Gaussian fits were only used to remove contaminat- 
ing emission (which represented only a fraction of the fiux 
in the main line) with non-parametric methods used to de- 
termine the line equivalent width, central wavelength and 
velocity widths. For the purposes of this fitting process the 
narrow emission lines and the broad Balmer lines (H(5, H7 
and H/3) were modeled as single Gaussians. When possible 
the number of independent parameters in the model was 
reduced by linking some together. For example, since the 
[O III] AA5007,4959 and narrow H/? emission arises from the 
same region of the QSO (the narrow line region, NLR) it 
seems reasonable to assume that the emitting gas will have 
similar velocity shifts and dispersions. 

The broad ultra-violet lines, i.e. from Mg II A2798 blue- 
ward, display emission line profiles with very broad bases 
which cannot be adequately modeled by a single Gaus- 
sian (Figure 1). They were therefore fitted with two com- 
ponents; a very broad Gaussian (full width at half maxi- 
mum, FWHM~ 10000 kms"^) and a narrower component 
(FWHM - 2000 - 4000 kms"^). These components were 
constrained to have the same central wavelength, with the 
exception of Lya as absorption to the blue side of the line re- 
sults in an asymmetric profile and it was therefore necessary 
to allow a velocity shift between the two components to fit 
the line profile. In all cases the best fit to the spectral feature 
was found using minimization techniques. See Groom et 
al. (2002) for further details on the line fitting procedure. 

No attempt was made to model and remove the broad 
Fc II emission underlying the Mg II feature. Our normali- 
sation procedure treats much of this feature as continuum 
emission, since there is no true continuum region close to 
the Mg II line, and it is therefore not possible to model and 
remove it from the normalised spectrum. The small Fc II fea- 
ture at ~ 29OOA just to the red of Mg II line, was modeled 
and removed. It was also not possible to deblend the Fe II 
emission lines underlying the H/3 spectral feature, the O IV] 



A1402 emission lines from the Si IV A1398 emission line and 
the Si III]] A1892 emission from the C III] A1909 feature (al- 
though the Al III A1857 emission was modeled and removed 
from C III]). 

Finally, the broad Balmer emission line H7 proved diffi- 
cult to deblend as it is contaminated by emission from both 
[Fe II] A4358 and [O III] A4364 as well as narrow H7 emis- 
sion. Since the [Fe II] and [O III] emission are within 6A 
of each other they are not resolved in the 2dF spectra and 
were therefore modeled as single narrow component centered 
between the two lines. The fit was further constrained by fix- 
ing the velocity width of the narrow H7 and the combined 
[Fc II] and [O III] Unes to that obtained for the [O III] A5007 
emission. 

Once we had modeled a spectral feature, the fits to the 
contaminating line emission were subtracted, leaving only 
the line of interest. Non-parametric methods (rather than 

the Gaussian fit) were used to measure the flux, velocity 
widths and central wavelengths of the emission line. These 
measurements were obtained for each feature both before 
and after the contaminating emission was subtracted in or- 
der to check for consistency. 

3.2 Measurement of line widths 

The emission lines widths were measured for the composites 
binned by magnitude and redshift (76 spectra in total; M—z 
composites hereinafter). Line widths have been measured 
using several different methods. These methods can be di- 
vided into two broad categories. The most common method 
(e.g. Boroson & Green 1992) is based on the full width of 
the line at some fraction of its maximum flux, e.g. the full 
width at half maximum (FWHM) or full width at zero in- 
tensity (FWZI), whereas the second method measures the 
velocity width outside which some fraction of the line flux 
falls, e.g.inter-percentile velocity widths (IPVs, e.g. Whit- 
tle 1985; Stirpe 1991; Stirpe, Robinson & Axon 1999). The 
most obvious difference between these two methods is that 
in the former method the emphasis is placed on the velocity 
width of the gas contributing to the core of the line emis- 
sion whereas the latter method takes into account the total 
flux in the line. The relative merits of the two methods are 
discussed elsewhere (sec, e.g. Robinson 1995; Stirpe 1991) 
and it is not clear which gives the most consistent estimate 
of the emitting gas velocity. For this study we measured the 
emission line width of our sample using both methods and 
comparison of the velocity widths measured for the same 
line using the two sytems shows they give consistent results. 

3.2.1 FWHM 

The FWHM is defined as the width of the line at the position 
where the line flux falls to half its peak value. The peak 
flux was taken to be the maximum flux within Ic ± O.Tbfa, 
where Zc and fa are the central wavelength and full width at 
half maximum obtained from the Gaussian fit to the data. 
In the case of lines fitted with multiple Gaussians fa was 
taken to be that of the broadest component. Interpolation 
between adjacent lA bins was used to estimate the FWHM 
to a fraction of an Angstrom and in the few cases where 
the line flux fell below half of the maximum flux and then 
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Table 1. Table showing results of rank Spearman correlation tests between the measured line width and the luminosity of the source 
for the M — z composite spectra. The number of data points, N, and the rank Spearman statistics, p, are shown together with the 
probability of that value of p being obtained by chance. These statistics are shown for the FWHM and IPV methods of measuring the 
velocity widths of the emission line and for the data before and after contaminating emission was subtra<;ted. 



(a) FWHM Method 




Entire Feature 


Contamination subtra<;ted 


Line 


N 


P 


P 


N 


P 


P 


Lya 


19 


0.789 


5.84e-05 


18 


0.564 


1.47e-02 


SilV 


25 


0.124 


5.55C-01 


25 


0.124 


5.55e-01 


CIV 


34 


0.669 


1.50C-05 


34 


0.679 


1.03C-05 


Willi T^OllllJ 


48 


0.618 


2.88C-06 


48 


0.426 


2.55C-03 


Mgll 


41 


0.574 


2.82C-04 


41 


0.574 


8.65C-05 


[NeV] 


23 


0.375 


7.83C-02 


23 


0.375 


7.83C-02 


[OTTl 


31 


0.139 


4.56e-01 


31 


0.139 


4.56e-01 


[Nelll] 


24 


0.258 


2.23e-01 


24 


0.283 


1.81e-01 


H5 


14 


0.473 


8.80e-02 


14 


0.240 


4.09e-01 


H-v 
nj 


20 


0.606 


4.62e-03 


20 


0.641 


2.34e-03 




15 


0.775 


6.90C-04 


15 


0.836 


1.04C-04 


fOTTTl 


17 


0.642 


5.45e-03 


17 


0.642 


5.45e-03 


(b) IPV Method 




Entire Feature 


Contamination subtracted 


Line 


N 


P 


P 


N 


P 


P 


Lya 


18 


0.740 


4.47e-04 


17 


0.716 


1.24e-03 


SilV 


30 


-0.054 


7.77e-01 


30 


-0.054 


7.77e-01 


CIV 


34 


0.678 


1.05C-05 


34 


0.623 


8.42C-05 


CllI]+SilII] 


48 


0.601 


6.35C-06 


48 


0.350 


1.47e-02 


Mgll 


43 


0.490 


8.41C-04 


43 


0.408 


6.61e-03 


[NcV] 


24 


0.712 


9.37C-05 


24 


0.712 


9.47C-05 


[Oil] 


31 


0.167 


3.68C-01 


31 


0.167 


3.68C-01 


[Nelll] 


25 


0.477 


1.59e-02 


25 


0.481 


1.50e-02 


H<5 


23 


0.447 


3.26e-02 


23 


0.242 


2.66e-01 


H7 


23 


0.752 


3.51e-05 


23 


0.611 


1.97e-03 


H/3 


17 


0.400 


1.12e-01 


17 


0.652 


4.57e-03 


[GUI] 


15 


0.907 


3.06e-06 


16 


0.906 


1.36e-06 



rose back above it, the largest possible width was used. This 
method has the advantage that unhke a Gaussian fit to the 
emission Une, no assumptions are made about the shape 
of the line profile but in low S/N data it is possible for 
the width of the line to be overestimated. The error in the 
FWHM is estimated from the error in the flux at the half- 
maximum positions and converted to an error in wavelengths 
using the gradient of the line profile at the half-maximum 
position. If the fiux in the lino fell below zero at any point 
between Ic ± 0.75/g the velocity width of that line was not 
measured. 



than just the core of the emission line. The error in the 
line width obtained via this method was estimated from the 
errors in the total integrated flux and the integrated flux 
at the 25 per cent and 75 per cent positions. The error in 
flux was then converted to an error in angstroms using the 
local gradient of the cummulative flux at that position. Since 
this method relies on the integrated line flux (rather than 
the peak flux) it can be used to measure widths for spectra 
with a lower SNR and we therefore obtained more line width 
measurements using this method than for the median width. 



3.2.2 IPV 

Here we use the integrated flux of the emission line to obtain 
an estimate of the line width. The cummulative flux in the 
line was found as a function of wavelength, beginning at the 
far blue side of the line {Ic — 1.5 x fa) with a cummulative 
flux = and ending at the far red side of the line {Ic + 
1.5 X fa), where the cummulative flux is equal to the total 
integrated flux of the line (F). The wavelengths at which 
the cummulative flux = 0.25F and 0.75F were found and 
their difference provided an estimate of the line width. Like 
the FWHM method described above, this method does not 
rely on any assumptions about the line profile and has the 
added advantage that all the flux in the line is used rather 



4 RESULTS 

4.1 Dependence of line velocity widths on QSO 
luminosity 

We derived the Spearman rank correlation coefficients for 

the correlation of line width (as measured by each method) 
with the source luminosity for each line (Table la and b). 
The emission lines widths were corrected for the instrumen- 
tal resolution {5\=9A for A=5500A in the observed spec- 
trum) by assuming that both the emission line and the in- 
strumental resolution were Gaussians. The small number of 
measurements with widths narrower than the instrumen- 
tal resolution were excluded as were velocity measurements 
with fractional errors of more than 20 per cent. We mear- 
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Table 2. Table showing results of partial Spearman correlation tests between the measured line width, the redshift and the luminosity of 
the sources contributing to the M — z composite spectra, taicing into account the correlation between redshift and luminosity. We show 
the partial rank Spearman correlation, p, and the probability of that value of p being obtained by chance. The number of data points 
used in this analysis, N, is the same as that listed in Table 1. The statistics are shown for the FWHM and IPV methods of measuring 
the velocity widths of the emission line and for the data before and after contaminating emission was subtracted. 



(a) FWHM Method 




Entire Feature 




Contamination subtracted 




Pv—z 


p 


Pv-L 


^v—L 


Pv-z 


p 

^v — z 


Pv-L 


Pv-L 


Lya 


-0.071 


7.72O-01 


0.756 


1.79e-04 


0.092 


7.17C-01 


0.508 


3.12C-02 


SilV 


—0.534 


5.94e-03 


0.383 


5.91C-02 


-0.534 


5.94C-03 


0.383 


5.91C-02 


CIV 


-0.674 


1.23e-05 


0.835 


8.55C-10 


-0.665 


1.77C-05 


0.835 


8.13e-10 


CIlI]+SilII] 


-0.271 


6.25e-02 


0.630 


1.58e-06 


-0.360 


1.21C-02 


0.536 


8.72e-05 


Mgll 


-0.036 


8.24e-01 


0.375 


1.57C-02 


-0.036 


8.24C-01 


0.375 


1.57e-02 


[NeV] 


-0.116 


5.98e-01 


0.270 


2.12e-01 


-0.116 


5.98e-01 


0.270 


2.12e-01 


[Oil] 


-0.269 


1.43e-01 


0.300 


1.02e-01 
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sured the correlation of the line width both before and after 
the contaminating emission was subtracted in order to de- 
termine whether our fitting procedure could be introducing 
spurious correlations into the data. 

A priori we chose a significant correlation to be one 
with a formal probability P < 1 per cent of being due to 
chance. For the M — z spectra significant correlations of line 
width with luminosity were found for the broad emission 
lines C IV, C III], Mg II and the Balmer lines H7 and H/3 
(Table 1) . With the exception of the H/3 emission line, these 
correlations were independent of the method used to mea- 
sure the line width and whether the contaminating emission 
had been subtracted. For H/3 only the IPV of the whole 
feature (H/3 plus the two [O III] lines) failed to show a sig- 
nificant correlation between luminosity and line width, most 
likely due to the contribution of the [O III] lines to the red 
side of the line. 

A highly significant correlation (P < 0.1 per cent) was 
also found between the width of the combined LyaN V fea- 
ture and the source luminosity but not in the N V-subtracted 
Lya line. This probably reflects the difficulties of accurately 
fitting the components in the LyaN V feature rather than 
a lack of a correlation between the width of Lya and lumi- 
nosity. The blue side of Lya emission line is subject to a 
considerable amount of absorption and it is probable that 



Lya is not well modeled using a combination of symmetrical 
line profiles, such as a Gaussians. Any errors in the fit to the 
Lya component would also cause errors in the fit to the N V 
component and the N V-subtracted Lya is therefore proba- 
bly not a good representation of the true line. Additionally, 
any variation in the relative strength of the Lya and N V 
lines, for example if N V became relatively stronger with 
luminosity, would also alter the FWHM of the feature. 

Aside from Lya, the only broad lines which do not ex- 
hibit a significant correlation between line width and source 
luminosity in the M — z spectra are Si IV and US. HJ is a 
relatively weak emission line and it is difficult to measure its 
width accurately. On the other hand Si IV docs not appear 
to vary significantly at all, either in fiux or velocity (Groom 
et al. 2002). 

Of the narrow emission lines, [O III] is the only line 
which exhibits a significant correlation with luminosity (Ta^ 
ble 1). 



4.2 Dependence of line widths on redshift 

To investigate whether there is any evidence for evolution in 
the AGN line width (i.e. narrower or broader lines observed 
at higher redshifts), we derived the Spearman rank correla- 
tion coefficients for a correlation of line width with redshift. 
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Figure 2. Velocity (■!;)-luminosity (Lj) plots for tiie FWH]V[ of the contamination subtracted spectra. Velocity is in units of kms~^ and 
Lb in units of erg s~^. The best fits of the form log('y) = a + 61og(Lj,) are also shown for each line (solid line for data which showed a 
significant correlation between velocity and luminosity; dashed lines for data which showed no significant correlation). We also show the 



best fit derived from regressing luminosity upon velocity (log(Lj,) 
equal to the number or degrees of freedom. 
© 0000 RAS, lUNRAS 000, 000-000 



c + (ilog(ii);dotted line). The errors bars are scaled such that x is 
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Table 3. Table showing the best fits for the velocity width, v (km s correlation with luminosity, L;, (erg s ^). A linear fit of the 
form log('y) = a + 61og(Lb) was applied to the data. 
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Significant positive correlations between redshift and veloc- 
ity were found for some methods of measuring the velocity 
of the emission linos (e.g. the FWHM velocity of Mg II, H7, 
and H/3; the IPV of [Ne V]) but unlike the luminosity cor- 
relations, no line showed a significant correlation between 
redshift and velocity for every method. 

A rank Spearman correlation test of our data reveals 
highly significant correlations {P < 0.01; generally P < 
10~® ) between luminosity and redshift for all the emission 
lines except Si IV and Lya and it is entirely possible that 
the handful of significant correlations we find between red- 
shift and velocity are due to a combination of the redshift- 
luminosity (z — L) correlation and that seen between lumi- 
nosity and velocity {v — L correlation) . We therefore derived 
the partial Spearman correlation coefficients (see Groom et 
al. 2002 for a detailed discussion of this analysis method) 
for the data. In this analysis we tost the hypothesis that any 
observed correlation between velocity and redshift is due to 
the z — L correlation and the v — L correlation with the co- 
efficients and their associated probabilities shown in Table 
2. For comparison we repeated this analysis for the v — L 
correlation, testing the hypothesis that the v — L correla- 
tion is due to a, V — z correlation and the z — L correlation 
(again shown in Table 2). We find that while the velocity 
is clearly correlated most strongly with luminosity for the 
C IV, C III] and Mg II emission lines, there is some evidence 
for a weak anti-correlation between velocity and redshift, 
significant only for the C IV emission line. If there is a weak 
trend for the emission lines to become narrower as redshift 



increases, this could easily be masked by the strong posi- 
tive v — L and L — z correlations, particularly as the range 
of redshifts over which the different emission lines arc ob- 
served is limited. The possible dependence of line velocity 
with redshift is discussed further in Section 5.2. None of the 
narrow emission lines exhibit any evidence for a correlation 
between redshift and velocity. 

We note that neither the H/3 or H7 lines exhibit a signif- 
icant v — L correlation in the partial Spearman correlation 
analysis although both lines show strong positive correla- 
tions using the standard Spearman correlation. They are 
also the only lines which exhibit a weak positive v — z corre- 
lation. This may be because the they arc not observed over 
a sufficiently large redshift range for the effects of the z — L, 
v — L and v — z correlations to be disentangled using this 
analysis (H/3, for example, is only observed in three redshift 
bins) . In any case it is clear from Figure 2 and the regression 
line analysis (Section 4.3) that both line widths of both H/3 
and H7 are correlated positively with luminosity. 

The interpretation of the partial Spearman correlation 
analysis for the Mg II emission line is also complex. Our re- 
sults indicate the Mg II line exhibits a weak correlation be- 
tween the FWHM velocity width and luminosity and a much 
stronger correlation between the IPV width and luminosity. 
The IPV widths also show a significant correlation with red- 
shift. We have been unable to correct for the Fe II emission 
underlying the Mg II line and if the relative strength of the 
Fe II features is a function of luminosity or redshift, this 
could cause the behaviour seen below in the Mg II line. For 
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example, an increase in QSO metallicity with redshift would 
increase the strength of the Fe II emission relative to that 
of the Mg II line and hence increase the width of the broad 
base (see Figure 1 and Section 3.1). 

4.3 Fits to the velocity-luminosity relationship 

For our composite spectra, the 6j luminosity Lb was esti- 
mated from the mean absolute bj magnitude of the spectra 
contributing to the emission line in each composite. The best 
fit to the relationship log(u) — a + blog{Lt) was obtained 
from a standard least-squares fit to the velocity width of 
the emission lines and QSO luminosity dependence. The pa- 
rameters a and b fitted to each line are given in Table 3 and 
a representative data-set is shown in Figure 2 (the FWHM 
velocity width of the contamination subtracted data). It is 
likely that the errors in the velocity widths have been un- 
derestimated in our analysis due to the method used to nor- 
malize the spectra, underlying Fe II features errors in the 
subtraction of the contaminating emission. The estimated 
errors on the emission line widths do not take these factors 
into account and hence the errors on a and b derived us- 
ing the standard method are probably an underestimate of 
the true error on the fit. We therefore used from the fit 
to rescale the errors on the velocity so that is equal to 
the number of degrees of freedom and hence achieve a more 
realistic estimate of the errors on a and b. 

As a further check of the correlation analysis performed 
above, wo refit the data, this time fitting a line of the form 
log(Lb) = c + dlog(v). For a perfect correlation between 
velocity width and luminosity (and no errors in either pa- 
rameter) the best fit lines of log(Lb) upon log(t;) will coin- 
cide with those of log(t;) upon log (Lb). For no correlation 
between v and L5 the lines of best fit will be perpendicular. 

The FWHM of the broad hydrogen linos, H/3, H7 and 
Lya show a slightly steeper dependence on luminosity than 
the other emission lines (Table 3). The best fit slopes to 
Mg II, C IV and Si IV lie in the range 0.04 < 6 < 0.1 whereas 
the slope of the hydrogen lines is in the range 0.1 < b < 0.2. 

Previous investigations into a possible correlation be- 
tween velocity width and QSO luminosity have produced 
ambiguous results. From observations of the H/3 emission 
line in individual sources KOO report an anticorrelation be- 
tween velocity and luminosity with a slope (in log-log space) 
of -0.27. On the other hand, several authors (e.g. Wandel 
& Yahil 1985; Joly et al. 1985; Boroson & Green 1992; 
Stirpe, Robinson & Axon 1999) have found the width of 
the Balmer lines (either FWHM, full width at zero intensity 
or full widths at some percentile of the maximum flux) to be 
weakly correlated luminosity. These studies were all based 
on samples containing <150 sources and it is possible that 
the inconclusive results may be duo to the intrinsic scatter 
in the v — L relationship as well as selection effects. 



5 DISCUSSION 

5.1 Broad emission lines and black hole mass 
estimates 

In reverberation mapping studies (e.g. Wandel et al. 1999, 
Onken & Peterson 2002) the rms spectrum of the H/? emis- 



sion line is often used to estimate the velocity dispersion of 
the gas. The rms spectra, formed from the rms difference 
between the individual spectra obtained during the moni- 
toring programme and the mean spectrum, measures only 
the components which vary and, since the emission from 
the NLR does not vary significantly on short time scales it 
should not be present in the rms spectrum. At first sight 
it therefore seems reasonable that the velocity dispersion of 
the BLR gas would be best represented by the H/9 veloc- 
ity width measured from the line after subtraction of the 
narrow emission line component. However, KOO found that 
there was little difference between the H/3 line width ob- 
tained from the rms spectrum and that measured from the 
mean spectrum without correcting for the narrow line emis- 
sion. Vestergaard (2002) also found that the width of the 
H/3 line in the rms spectrum was more similar to that of 
a single-epoch observation before any correction for NLR 
emission. This implies that much of the narrow core of the 
H/3 emission line (Figure 1) is from gas associated with the 
BLR (rather than the NLR) and wc therefore use the H/3 line 
width from the composite spectra before correction for any 
narrow line emission in the analysis which follows. For the 
non-Balmcr lines, wo use the FWHM line widths measured 
after contaminating emission has been subtracted. 

A secondary issue is the underlying Fe II emission. As 
discussed in Section 3.1, many of the Fe II emission features 
have been treated as continuum, particularly the large fea- 
ture underlying Mg II, and it is therefore not possible for 
us to correct the emission lines for the underlying Fe II. 
This may impact on our measurement of the Mg II FWHM. 
In their study, McLure & Jarvis (2002) used Mg II FWHM 
measurements from spectra taken with the International Ul- 
traviolet Explorer (lUE) which had been corrected for the 
underlying Fe II emission. They found a linear almost one- 
to-one, correlation between the FWHM measured for the H/3 
and Mg II lines whereas wc find that the Mg II emission is 
generally broader than the H/3 emission. On the other hand, 
some of the H/3 FWHM measurements used by McLure & 
Jarvis (2002) have been corrected for narrow H/3 emission 
and in these cases it is reasonable to asume that their H/3 
FWHM would be broader than ours. 

There is also an Fe II emission line to the red of H/3 
(~ 4930 A, see e.g. Vanden Berk et al. 2002) which we have 
not removed before measuring the H/3 FWHM. Vestergaard 
(2002) found that it was not necessary to correct for the 
Fell emission unless the Fell emission from the object was 
unusually strong, and as there is no evidence for strong Fe II 
emission in our composite spectra, we do not believe this 
feature will have a significant impact on the H/3 FWHM. In 
any case, we note that KOO did not correct for Fe II emission 
when measuring the H/3 FWHM from their mean spectra. 

Having determined which measurements of the line 
widths most accurately represent the velocity dispersion of 
the gas, wo need to obtain an estimate of the radius of the 
BLR. From reverberation mapping experiments, KOO found 
the BLR radius, r oc Lsioo where isioo is the continuum lu- 
minosity at 5IOOA (defined as AL5100 in KOO and = 3.56Lb 
for our data) . This result depends on the cosmology assumed 
and the analytic method used to determine the fit. KOO 
assumed qo = 0.5, Ao = and Ho = 75 kms~^ Mpc~^ 
and used the fitexy linear regression method (Press et al. 
1992), which takes into account errors in both L and r, 



© 0000 RAS, MNRAS 000, 000-000 



10 E. A. Corbett et al. 




45 45.5 46 46.5 44 45 46 43.5 44 44.5 

Log(L) Log(L) Log(L) 



Figure 3. The mass (M)-luminosity (L5100) relationship for the different broad emission Hnes. Mass is given in units of Mq and 
luminosity in erg s~^. Different colours reflect the different redshift ranges as shown in the key for each plot and the error bars shown 
have been scaled such that is equal to the number of degrees of freedom. The YonX (solid black line) and BCES (green dot-dash line) 
fits to the Mbh — L relationship for each line are also shown. 



to fit their data. Re-analysing the KOO data, assuming the 
cosmology adopted here {Qo = 0.3, Ao ~ 0.7 and Ho — 
70 kms~^ Mpc~^), Netzer (2003) finds a similar relation- 
ship with 
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Using the BCES bisector method (Akritas & Bershady 
1996) to fit the same data, Netzer (2003) obtains a somewhat 
flatter relationship, r oc isioo*"'^^, although the two fits are 
formally consistent. Similar results are obtained by McLure 
& Jarvis (2002) who find r oc L°foo^°-^° for our cosmology 
and the BCES bisector method. Here we have adopted r oc 

^5100 (Equation 2) but we note that a flatter r I/sioo 

relationship will yield smaller black hole masses at higher 
luminosities. 

The black hole mass is then given by (KOO): 
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'light days" 1000km s-i ' """^ 

where v is the FWHM of the H/3 emission line. The 
BLR clouds are assumed to be in random orbits about the 
black hole and thus in Equation 3 the velocity dispersion is 
assumed to be (\/3/2)« (Netzer 1990). If the orbits of the 
BLR clouds are not random but confined to disks, as sug- 
gested by e.g.. Young et al. 1999, McLure & Dunlop (2002), 



the velocity dispersion is {3/2)v, and the black hole masses 
are a factor of 3 larger than we derive. 

The black hole masses for selected lines were estimated 
using Equations 2 and 3 and are shown in Figure 3 plot- 
ted against Lsioo- As expected, a strong correlation can be 
seen between BH mass (M) and the luminosity, L5100, for 
each line with remarkably little scatter, generally less than 
0.3 dex. The best fit of the form log(M) = a + /31og(L) 
for each line is shown in Figures 3 and 4 with the parame- 
ters for each fit given in Table 4. The data were fitted using 
both the standard linear regression method (YonX) and the 
B CSS method. From Table 4 it is clear that the two methods 
agree to within the errors (see also Figure 3). For the lin- 
ear regression method the errors in the fit parameters were 
evaluated as described in the previous section, using the 
of the fit to derive the errors in the mass estimate and hence 
a and /3. We note that the linear relationship between Mbh 
and L5100 derived for the H/3 line is similar to that obtained 
for the Lya emission line in this analysis. Fits to the other 
emission lines are slightly offset from, and fiatter than, the 
fit to the hydrogen lines (Figure 4). Nevertheless, for a given 
continuum luminosity, the black hole masses obtained from 
the different emission lines typically lie within ±0.3 dex of 
each other. 

Reverberation mapping experiments (e.g., Peterson & 
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Table 4. The best fits to the relationship between the black hole mass, M, with the monochromatic luminosity, Z/510O) in units of erg 
s^^ . Black hole mass is in units of Mq. A linear fit of the form log{M) = a + l3log{L5ioo) was applied to the data, where f} is the slope 
of the relationship and a is a constant. S is the normalisation constant applied to the data to force the weighted mean of the line to lie 
on the H/3 YonX relationship and has errors of ~ O.ldcx 
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Figure 4. The mass (M)-luminosity (L5100) relationship for se- 
lected broad emission lines, each line shown in a different colour. 
The best BCES fit to the Mbh — L relationship for each line is 
also shown. Error bars are not plotted for clarity. 



Wandel 1999, 2000; Onken & Peterson 2002) find that diflter- 
ent emission lines respond to changes in the UV continuum 
with different time lags. This indicates that the BLR is strat- 
ified, with the broader emission lines (such as C IV) origi- 
nating closer to the central engine than the Balmer lines. 



The empirical r — I/5100 relationship (Equation 2) was ob- 
tained by KOO from the H/3 and Ha emission lines and thus 
its use to determine black hole masses via Equation 3 will 
lead to an overestimate of the black hole mass for C IV. To 
obtain more reliable estimates of M for the UV lines it is 
necessary to introduce a "calibration" factor to equation 2 
to correct for the difference between the radius of the gas 
emitting the broad UV lines, such as C IV, and that emitting 
the H/3 lines. This has already been attempted by Vester- 
gaard (2002) who compared the BH masses estimated from 
reverberation mapping of the H/3 emission line with those es- 
timated from the C IV FWHM in the same object assuming 
that, for the C IV line, 

il350 



FWHMc iv.2f 
- 1000km s-i > 4044erg s-i ^ 



(4) 



where I/1350 is the continuum luminosity at 1350A (i.e. 
ALi3so in Vestergaard 2002) and a is a constant. She found 
a = io6-2±0-45 which implies that the radius of the C IV emis- 
sion, rc IV ~ 0.5'^Q'\rHi3 (assuming that the QSO spectra 
have a spectral index of —0.5) although the errors in a, and 
hence this result, are large. 

To determine the calibration factor for our data we 
re-normalised the mass-luminosity relationship so that the 
weighted mean of the mass derived from each line lies on the 
H/3 (YonX) relation, noting that an extrapolation of the H/3 
mass-luminosity relation to higher luminosities is entirely 
consistent with that observed for Lya. However, because 
our estimates of black hole mass are based on the under- 
lying assumption that the KOO radius-luminosity relation 
r oc l'^-^^ holds for all emission lines, we chose not to ad- 
just the slopes of the observed mass-luminosity relation for 
each emission line. The mass-luminosity relationships were 
therefore normalised by the addition of a constant offset, S 
(Table 4). Assuming that this offset is due to the emission 
lines being emitted from gas situated at different radii, 5 can 
be converted to provide an estimate of the radius at which 
the UV lines are emitted relative to that of the H/3 emission 
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Figure 5. The mass (M)-luminosity (L5100) relationship for the 
different lines normalised to the H/3 relationship. Each colour cor- 
responds to a different line. The best BCES fit to the data as a 
whole is shown (solid line) as well as that to the If/3 and C IV lines 
alone (dotted line). Error bars are not plotted in the interests of 
clarity. 

(radius = lO'^rn/j). We find that the rc iv = 0.6lg irH,3 
(Table 4) which is consistent with that found by Vester- 
gaard (2002). From their reverberation mapping study of 
NGC 3783, Onken & Peterson (2002) also find that the C IV 
emission line lag is approximately half that of the H/3 line 
although the errors on this measurement are also large. 

The calibration factor S does not address the variation 
in the slopes of the mass-luminosity relationships for the 
different lines, which, if significant, must result from a dif- 
ference in the measured line width-luminosity relation for 
different species. Since it is difficult to see how different 
emission lines could give rise to a different mass-luminosity 
relations (a line-independent quantity) based on the same 
set of QSO composite spectra, the most straightforward in- 
terpretation of this observation is that the different slopes 
are caused by the contaminating emission blended with the 
lines. As we have already discussed in this and previous sec- 
tions, the Mg II and C III] emission lines are blended with 
emission from other elements (Fe II and Si III] respectively) 
which would effectively broaden the measured FWHM of the 
lines, making our calculation of S unreliable. The blended 
emission could also alter the slope of the log(t;) — \og{L) re- 
lation, particularly if the relative contribution of the lines 
vary with luminosity or redshift. We note that the slope 
of the log(M) — log(L) relationship for H/3 and C IV are 
statistically identical given their errors (Table 4) and that 
these lines also suffer the least amount of contamination. 
We cannot, however, rule out the possibility that the radius- 
luminosity relation may differ slightly from line to line. 

The re-normalised mass-luminosity relations for all lines 



are plotted in Figure 5 and we fit an overall mass-lumirrosity 
relation of the form: 

log(AfBH) = 0.93(±0.05) log(L5ioo) - 33.5(±2) (5) 

to all emission lines (solid line. Fig 5). However, mindful the 
problems associated with contaminating emission, we also 
fit a mass-luminosity relation to the H/3 and C IV emission 
lines only (Fig. 5; dotted line) since these lines are the least 
affected by contaminating emission and should therefore give 
the most reliable estimate. For the H/3 and C IV lines alone, 
we find 

log(i\/BH) = 1.00(±0.1) log(L5ioo) - 36.8(±4.5) (6) 

which is slightly steeper than that found for all the emission 
lines combined. 

The errors on the Mbh — L relations are larger than 
those quoted in Table 4. In estimating the errors on these 
relations we have taken into account the errors on the ini- 
tial fit to H/3 which, since it is used to normalise the mass 
estimates from the other emission lines, contributes to the 
overall errors on the Mbh — L relation. The error in the fit 
to the H/3 and C IV lines is larger than that obtained for 
all the emission lines as we have fewer data points. Note the 
errors on the gradient and the intercept in Equation 5 and 6 
are not independent and have a correlation of nearly -1. The 
relations we obtain are in reasonable agreement with the sta- 
tistical analysis of Netzer (2003), who finds Mbh oc L^f^"-^^ 
for tblr oc 1/1350 (where Z/1350 is the luminosity at 1350A) 
for a sample of over 700 objects (Table 1; Netzer 2003). 

An important caveat on the Mbh — L relationship given 
in equations 5 and 6 is that throughout our analysis we have 
used the relationship tblr oc Lgj^oo discovered by KOO and 
later refined by McLure & Jarvis (2002) and Netzer (2003). 
As discussed previously, the value of 7 is uncertain and de- 
pends on the statistical method used to fit the data, with 
the fitexy method and the BCES method giving values of 
0.68 ± 0.03 and 0.58 ± 0.12 respectively for the same data. 
Although there is a correlation between velocity and lumi- 
nosity (Section 4) the slope of the log(i') = a + 61og(L5ioo) 
relationship is quite flat for most of the emission lines (e.g, 
b = 0.17 for H/3, 0.1 for C IV, 0.05 for Mg II; Table 3) and 
the slope of the log(MBH) — log(L5ioo) relation is therefore 
dominated by the value assumed for 7 which is currently 
known to an accuracy of only ±0.15. For example, if we 
adopt 7 — 0.58 in Equation 2, we find Mbh oc L'^foo for 8-11 
lines and Mbh oc Lsioo for the H/3 and C IV lines alone, 
again consistent with Netzer (2003). 

Additionally, and perhaps more importantly, it is not 
known whether 7 is constant across all luminosity ranges. 
The sources in KOO's sample all have I/5100 < lO**" erg s~^, 
with only 4 (out of a total sample of 34) having L5100 > lO**^ 
erg s~^. Our mass estimates for the UV lines are therefore 
based on the assumption that the tblr oc I/5100 relation can 
be safely extrapolated to higher luminosities. Netzer (2003) 
found that eliminating the 3 objects with L5100 < lO'*'^ erg 
s"^ from KOO's sample resulted in a -^{BCES) = 0.71 ±0.21 
and '■/{fitexy) = 0.69 ± 0.03 whereas eliminating the seven 
objects with L5100 < 10*^ '^ erg s"^ resuhed in ^{BCES) = 
0.58±0.19 and y{fitexy) = 0.74±0.04. Reverberation map- 
ping of multiple line species in the same object and extend- 
ing the luminosity range of reverberation mapped sources 
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would clearly be of great benefit in eliminating some of un- 
certainties with this analysis. 

5.2 Narrow emission lines and black hole masses 

A number of studios have found that the stellar velocity dis- 
persion of the host galaxy bulge is related to the bulge lumi- 
nosity which is, in turn, related to the mass of the black hole 
(e.g. Magorrian et al. 1998; Ferraxese et al. 2001; Tremaine et 
al. 2002). Recently, Shields et al. (2002) and Boroson (2003) 
investigated the possibility of using the velocity width of 
[O III] as a substitute for the stellar velocity dispersion when 
determining the black hole mass. They both found that the 
velocity width of [O III] is correlated with black hole mass, 
assuming that the black hole mass oc U'v^ where L is the 
QSO luminosity at 5100A, v is the velocity width of the H/3 
emission line and 7 = 0.5 and 0.7 for Shields et al. (2002) 
and Boroson (2003) respectively. We find that the velocity 
width of [O III] is correlated with the luminosity of the QSO 
(Table 1) with v oc L° i~° is (Table 3; FWHM values as 
the IPV values arc probably contaminated by Fe II omis- 
sion either side of the [O III] lines), similar to that found 
for H/3. This raises the intriguing possibility first, that the 
NLR in QSOs is sufficiently close to the black hole for its 
motion to be dominated by the gravitational potential of 
the black hole, and secondly, that the narrow emission lines 
may be used to derive the black hole mass. However, none of 
the other narrow emission lines in our study ([O II], [Ne V] 
and [Ne III]) display the same strength of correlation be- 
tween line width and QSO luminosity as [O III] and given 
the anti-correlation between narrow lino strength and con- 
tinuum luminosity (Groom et al. 2002) the use of narrow 
lines in estimating the black hole mass may be limited to 
low luminosity and hence low mass objects. 

5.3 Evolution of the black hole 
mciss-to-luminosity ratio 

Figure 3 shows the mass-luminosity relation for the indi- 
vidual broad lines separated into redshift bins and appears 
to demonstrate a lack of evolution in the relation between 

black hole mass and luminosity with redshift. However, the 
C IV line does show a significant anti-correlation between 
velocity and redshift in our partial correlation analysis (Ta- 
ble 2) which may be evidence of a weak redshift evolution of 
line width. No other lines show significant correlations with 
redshift. 

In order to make a quantitative estimate of the amount 
of redshift evolution we fit a simple power law model to the 
L — z data points for each emission line. Wo first assume 
that black hole mass is related to luminosity and redshift by 
independent power laws such that 

Mbh(.?)ocL(«)"(1 + «)'^, (7) 

and that secondly, the radius-luminosity relation is also de- 
scribed by a power law: 

R(xL{zy. (8) 

We assume that the radius-luminosity relation (Equa- 
tion 8) does not to evolve with redshift as it should be pri- 
marily driven by the physics of photoionization, which we 
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Figure 6. The evolution of C IV velocity width. We plot the 
re-normalized (to a fixed absolute magntude of Mb = —25.0) 
velocity vs. (I -|- z) in log-log space. The error bars have been re- 
scaled such that is equal to the number of degrees of freedom. 
A decline in velocity can be seen at the highest redshifts. 

would not expect to change with redshift (although this may 
not be the case if QSO metallicity was found to evolve). As- 
suming Keplerian motion, we have 

\og{v) = log(L) + ^ log(l + z)+ Const. (9) 

This relation is then fitted to the observed velocity- 
luminosity data for each emission line individually (as we 
do not wish to assume that the radius-luminosity relation is 
the same for each emission line). To obtain a more realistic 
estimate of the errors on the fitted parameters given a large 
intrinsic scatter in the data, we re-scale the errors on indi- 
vidual points such that is equal to the number of degrees 
of freedom. The results of this model fitting are shown in 
Table 5. The only line which shows clear evidence (> 3a) 
for evolution is C IV. However, Si IV and C III] also show 
weak evolution in the same sense. All of these lines have best 
fit values of /3 which are negative, implying that for a given 
luminosity, the measured velocities are lower at higher red- 
shift. By contrast, the value of f3 for the Mg II line is entirely 
consistent with no evolution at all and the Balmer lines (H/3 
and H7) all show positive values of (3, although these are 
not significant. We do not attempt to model the redshift 
evolution of the Lya line, as the contribution of Lya for- 
est absorption and N V emission makes any interpretation 
highly uncertain. To show more clearly the evolution in the 
C IV velocity width, we plot log(v) vs. log(l-|-«) for this line, 
first renormalizing the velocity points to a constant absolute 
magnitude (based on the above fits) of Mb ~ —25.0 (Fig. 
6). The error bars in Fig. 6 have been rescaled such that 
is equal to the number of degrees of freedom. A weak, but 
significant, decline in velocity can be seen, particularly at 
log(l -I- 2) > 0.5. 

The evolution seen in the higher redshift lines is in the 
sense of higher redshift QSOs having smaller velocities for 
the same luminosities, which might be expected if QSO lu- 
minosity evolution is driven by a decline in fueling rate at 
lower redshift. However, the amount of evolution is small 
compared to the strong luminosity evolution shown by QSOs 
(typically ~ (1 + .z)'^). At « ~ 2 - 3 the QSO evolution 
slows down and appears to turn over. Assuming pure lumi- 
nosity evolution (PLE), Groom et al. (2003) finds that for 
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the full 2QZ catalogue the evolution of L* , the characteris- 
tic luminosity of the QSO population, is well described by 
L*{z) = L*(0)10'=i^+*==^^' with fci = 1.17 and fea = -0.21 
at redshifts less that z ~ 2.3. This model implies that L* 
reaches an extrapolated maximum at 2; ~ 2.8. If black hole 
masses do not evolve, and the luminosity evolution is purely 
due to other influences (e.g. fueling rate), then at a given lu- 
minosity Mbh(2) = A^bhCO). Thus, using the PLE fit from 
Boyle et al (2003) and assuming that M oc we find 

Mbh{z) _ Mbh(O)^^ 

L"{z) ~ L«(0) ■ ^ ' 

In Figure 7 wo compare this model to the evolutionary slopes 
found in our fits. Our PLE model (solid line) is steep at low 
redshift, but flattens ofi^ at 2 ~ 2 — 3. Here we assume that 
a = 0.93, as found in Section 5.1 for all broad emission lines, 
but the shape of the model curve changes little if a = 1.0, as 
found for the H/3 and C IV lines alone (Section 5.1). We plot, 
as shaded regions, the range of allowable (±lcr) slopes for 
each of flve lines (H/3, Mg II, C III], C IV and Si IV), over the 
redshift ranges that they are measured. Each is normalized 
to the M{z)/L"{z) ratio at the mean redshift of the line. At 
low redshift, the H/3 line has a positive slope which is con- 
sistent with /3 = 0, but not with the strong evolution of the 
model (solid line). At higher redshift the Mg II and C III] 
linos arc close to f3 — and although this is inconsistent 
with our model, we note that both these lines are contam- 
inated by blended emission which could be distorting their 
behaviour. Only the highest redshift lines, C IV and Si IV 
show evolution which is comparible to that shown by the 
model, partly as they have higher values of /3 and partly as 
the model starts to turn over at high redshift. 

It appears that only aX z > 2 does the measured evolu- 
tion in the luminosity- velocity relation compare to the above 
model. Over the majority of the redshift range covered in 
this investigation evolution in the luminosity-velocity rela- 
tion is insigniflcant. This conclusion has significant impli- 
cations for QSO evolution, implying little evolution in the 
mass-luminosity relation for super-massive black holes over a 
wide range in redshift. Assuming no evolution in the mass- 
luminosity relation wo can transform the QSO luminosity 
function (LP) and its evolution with redshift into a mass 
function (MP) . In Pig. 8 we show the evolution of the QSO 
mass function based on the LP derived from the 2QZ-I-6QZ 
samples (Groom et al. 2003). We immediately note that the 
luminosity (mass) evolution of the QSO LP cannot bo ex- 
plained by decline in the fueling rate - either in continu- 
ous model (long-lived QSOs) or in a duty cycle (short-lived 
QSOs), since this would imply that a QSO at a given rela- 
tive position on the LP (MP) evolving from high black hole 
meiss at high redshift to low black hole mass at low redshift. 

Rather, the evolution of the black hole mass function 
suggests that QSOs can only be active once, possibly associ- 
ated with formation of the black hole itself during spheroid 
formation. Under this interpretation, the evolution of the 
LP is driven by the efficiency of black hole formation which 
drops rapidly towards low redshift. This picture is qualita- 
tively consistent with the recent theoretical work of Kauff- 
mann & Haehnelt (2000), although it places tight constraints 
on the fraction of QSOs which exhibit a duty cycle of activ- 
ity in such models. 

It also implies that the Gebhardt et al. (2002) relation 
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Figure 7. Tlie best fit evolutionary power law relationships. Tiie 
evolution predicted by the luminosity evolution of Boyle et al. 
(solid line) is compared to the slopes found in our fits (shaded 
regions). The horizontal width of the individual shaded regions 
represents the redshift range covered by our data for each emission 
line and the height of the regions represents the la error on (3 
(Table 5). The vertical positioning of the shaded regions is such 
that they intersect the model at their mean redshift (see text for 
a more detailed discussion). 

between black hole mass and galaxy velocity dispersion will 
evolve toward higher ratios at higher redshifts. The black 
hole mass will be 'frozen' when black hole formation and 
subsequent QSO activity ceases while the mass/velocity dis- 
persion of the host galaxies will increase due to hierarchical 
merging to give the relationship seen at the present day. 

At 2 ~ 2 and greater, there does seem to be some evolu- 
tion of the velocity-luminosity relation, as shown by the C IV 
(Pig. 6) and Si IV lines. It appears that there is a general 
decline in velocity for a given luminosity at log(l -\- z) > 0.5 
{z > 2.2). The measured rate of evolution is no difi'erent 
if we remove the highest redshift bin (2.75 < z < 3.0) 
and re-fit the data. The fact that the strongest evolution in 
the velocity-luminosity relationship is found in the redshift 
regime corresponding to the peak of QSO activity suggests 
that as L* declines at high redshift [z > 2.5), the black hole 
mass-luminosity relation may be changing more rapidly than 
at lower redshift. 



6 CONCLUSIONS 

We have used composite spectra generated from more than 
22000 QSOs observed in the course of the 2dP and 6dP QSO 
redshift surveys to investigate the relationship between the 
profiles of emission lines and luminosity. We find that the ve- 
locity width of the broad emission lines H/3, H7, Mg II, C III] 
and C IV show a positive correlation with the contiimum lu- 
minosity (as measured in the bj filter) but the slope of this 
dependence is much flatter for the UV emission lines (Mg II, 
C III] and C IV ) than for the hydrogen lines. This correlar 
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Table 5. The best fit evolutionary power law slopes. 
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Figure 8. The Black Hole Mass Function for QSOs and its evolution with redshift. 



tion is statistically significant with a confidence level > 99 
per cent and is independent of tlie metliod used to measure 
the emission line widths. From partial Spearman correlation 
analysis we find only weak evidence for an anti-correlation 
between redshift and velocity. Of the narrow emission lines, 
[O II], [Ne III], [Ne V] and [O III], only the [O III] A5007 line 
exhibits a correlation between line width and luminosity. 

Assuming that the gas in the BLR is in near-Keplerian 
orbits and that the radius of the BLR, r oc L°-^^ (KOO;Netzer 
2003), where L is the continuum monochromatic luminos- 
ity at 5IOOA, we have used our measurements of the emis- 
sion line velocity widths to derive estimates of the average 
black hole mass in the composite spectra. We find that for 
a given luminosity the black hole masses obtained from the 
different the emission lines agree to within 0.6dex. However 
the scatter in the mass-luminosity relationships for the in- 
dividual lines is much smaller (< 0.3 dex) and given that 
the different emission lines may be emitted from different 
radii, producing an offset between the masses estimated from 
different emission lines, we use the H/3 Mbh — L relation- 
ship to calibrate the other emission lines. We find an over- 
all black hole mass - luminosity relationship of the form 
Mbh oc . It is possible that the FWHM measure- 

ments of several of the emission lines are overestimated as 
they are blended with weaker lines from different elements 



which are difficult to remove reliably. We have therefore 
also derived the mass-luminosity relationship for the H/3 
and C IV emission lines only, as these are the lines least ef- 
fected by contaminating emission, and find Mbh oc 7;,i o±o i, 
OurA/BH — L relationships depend heavily on the slope of 
the log(r) — log{L) relationship which is currently known to 
an accuracy of only ±0.15. Using the midpoint of the de- 
rived Mbh — L powerlaws and taking into account the error 
in the slope of the log(r) — \og{L) relationship, we conclude 
that Mbh oc L''•^^±0■l^ 

We find only weak evidence of evolution in the velocity- 
luminosity correlation, which places strong constraints on 
the evolution of black hole masses for a fixed QSO luminos- 
ity; M oc (I + z)^ , P < 1. This is much smaller than that 
required for the luminosity (mass) function of QSOs to be 
due to a single population of long-lived sources. 



ACKNOWLEDGMENTS 

The 2dF QSO Redshift Survey was based on observations 
made with the Anglo-Australian Telescope and the UK 
Schmidt Telescope. We warmly thank all the present and 
former staff of the Anglo- Australian Observatory for their 
work in building and operating the 2dF facility. 



© 0000 RAS, MNRAS 000, 000-000 



16 E. A. Corbett et al. 



REFERENCES 

Akritas M.G., Bershady M.A., 1996, ApJ, 470, 706 

Baldry I. K. ct al, 2002, ApJ, 569, 582 

Baldwin J. A., 1977, ApJ, 214, 679 

Blanford R.D., McKcc C. F. 1982, ApJ, 255, 419 

Boroson T. A., Green R. F., 1992, ApJ, 80, 109 

Boroson T. A., 2003, ApJ in press (astro-ph/0211372) 

Boyle B. J., Shanks T., Groom S. M., Smith R. J., Miller L., 

Loaring N., Hcymans G., 2000, MNRAS, 317, 1014 
Boyle B. J., 1990, MNRAS, 243, 231 

Boyle B. J., Shanks T., Peterson, B. A., 1988, MNRAS, 235, 935 
Gristiani S., Vio R., 1990, A&A, 227, 385 

Groom S. M., Smith R. J., Boyle B. J., Shanks T., Loaring N. S., 
Miller L., Lewis 1. J., 2001, MNRAS, 322, L29 

Groom S. M., Rhook K., Gorbett E. A., Boyle B. J., Netzer H., 
Loaring N. S., Miller L., Outram P. J., Shanks T., Smith R. J., 
2002, MNRAS, 337, 275 

Groom et al. 2003, in preparation 

Ferrarese L., Podge R. W., Peterson B. M., Merritt D., Wandel 

A., Joesph, G. L., 2001, ApJ, 555, L79 
Francis P. J., Hewctt P. G., Foltz G. B., Ghaffee F. H., Weymann 

R. J., Morris S. L., 1991, ApJ, 373, 465 
Gebhardt K., et al. 2000,ApJ,543, L5 

Green P. J., Forster K., Kuraszkiewicz J., 2001, ApJ, 556, 727 
Joly M., GoUin-Souffirin S., Masnou J.L., Nottale L., 1985, A&A, 
152, 282 

Kaspi S., Smith P. S., Maoz D., Netzer H., Jannuzi B. T., 1996, 
ApJ, 471,L75 

Kaspi S., Smith P. S., Netzer, H., Maoz, D., Januzzi B. T., Giveon 

U., 2000, ApJ, 533, 631 (KOO) 
KauflFman G., Haehnelt G., 2000, MNRAS, 311, 576 
KroUck J., ApJ, 551, 72 
Magorrian J., et al. 1998, AJ, 115, 2285 
Maoz, D. 2002 astro-ph/0207295 

Miyoshi M., Moran J., Hernstcin J., Grecnhill L., NaJsai N., Dia- 
mond P., Inoue M., 1995 Nature, 373, 127 
McLure, R. J., Dunlop, 2002, MNRAS, in press 
McLure, R. J., Jarvis, M. J., 2002, MNRAS, in press 
Netzer H., 1990 in Active Galactic Nuclei, Saas Fee Ad- 
vanced Course 20, ed. T. J. -L. Gourvoisier & M. Major 
(Berlin: Springer), 57 
Netzer H., 2003, ApJ, 583 L5 
Onken C. A., Peterson B.M. 2002, ApJ, 572, 746 
Peterson B. M., Wandcl A., 1999, ApJ, 521, L95 
Peterson B. M., et al., 2000, ApJ, 542, 161 

Press W. H., Teukolsky S.A., VctterUng W.T., Flannery B.P., 
1992, Numerical recipes (Second ed,; Cambridge: Cambridge 
Univ. Press), 660 

Robinson, A., MNRAS, 272, 647 

Schlegel D.J., Finkbeiner D.P.,Davis M., 1998, ApJ, 500 525 
Smith J.E, Robinson A., Young S., Gorbett E.A., Axon D.A., 

Gianuzzi M.E., 2002, MNRAS, in press 
Stirpe G., 1991, A&A, 247, 3 

Stirpe G., Robinson A., Axon D. J, 1999, in Structure and Kine- 
matics of Quasar Broad Line Regions, ed. C. M. Gaskell, W. 
Brandt, M. Dietrich, D. Dultzin-Hacyan, Eracleous, 175, p. 
249 

Shields G. A., Gebhardt K., Salviander S., Wills B. J., Xie B., 
Brotherton M.S., Yuan J., Dietrich M., 2002, ApJ submitted 
(atsro-ph/0210050) 

Tanaka Y., et al. 1995, Nature, 375, 659 

Tremainc et al., 2002, ApJ, 574, 740 

Wandel, A., Peterson, B. M., Malkan, M. A., 1999, ApJ, 526, 579 
Wandel, A., A., 1985, ApJ, 295, LI 
Whittle,1985, MNRAS, 213, 1 

Wilms J., Reynolds G. S., Begelman M. C, Reeves J., Molendi 
S., Staubert R., Kendziorra E., 2001, MNRAS, 328 L27 



Woo J. H., Urry G. M., 2002, ApJ, 579, 530 

Young S., Gorbett E.A., Hough J.H., Gianuzzi M.E., Robinson 

A., 1999, MNRAS 
Vestergaard, M. 2002, ApJ, 571, 733 



© 0000 RAS, MNRAS 000, 000-000 



